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We study the dimer-dimer scattering length 04 for a two-component Fermi mixture in which the 
different fermions have different masses rrif and mj.. Fhis is made in the framework of the exact 
field theoretical method. In the large mass ratio domain the equations are simplified enough to lead 
to an analytical solution. In particular we link 04 to the fermion-dimer scattering length ag for the 
same fermions, and obtain the very simple relation 04 = as/2. Fhe result 04 ~ 03/2 is actually 
valid whatever the mass ratio with quite good precision. As a result we find an analytical expression 
providing a^ with a fairly good precision for any masses. Fo dominant orders for large mass ratio it 
agrees with the literature. We show that, in this large mass ratio domain, the dominant processes 
are the repeated dimer-dimer Born scatterings, considered earlier by Pieri and Strinati. We conclude 
that their approximation, of retaining only these processes, is a fairly good one whatever the mass 
ratio. 

I. INTRODUCTION 

Ultracold atoms are a remarkable playground for a number of other fields of physics, such as condensed matter 
physics, nuclear physics and astrophysics. This is due to the simplicity at very low temperature of the effective 
interaction between atoms and moreover to the experimental ability to choose almost at will the corresponding 
parameters jT]. The resulting physical systems correspond often to simple limiting situations of high interest in 
these other fields. The experimental realization of the BEC-BCS crossover in Fermi gases is a wonderful example 
of the flexibility available in cold gases. As a result of the ability to vary the effective interaction in a very wide 
range, it has been possible to go from the weakly attractive regime between fermions, where at low temperature a 
BCS-like condensate arises, to the strongly attractive regime where molecules, or dimers, form and these composite 
bosons undergo the standard Bose-Einstein condensation at low temperature. This crossover has provided an explicit 
demonstration of the deep-seated connection between these two kinds of superfiuids, which had been realized long 
ago [2HS]- At the same time it displays, at unitarity where dimers are just appearing, a somewhat new example of 
superfluidity. 

When the interaction strength is just beyond the one corresponding to unitarity so that dimers are just forming, the 
physical situation is quite complex since the dimer size is very large and they overlap very strongly in the gas leading 
to a complicated many-body problem, at higher temperature in the normal state as well as in the superfluid state 
found at lower temperature. However when the interaction strength is increased toward the strong coupling regime, 
the dimer size becomes very small and their overlap becomes negligible. In this case the composite nature of the 
dimers becomes irrelevant and physically one deals with a simple boson gas. Just as for ultracold Bose gas, the only 
relevant parameter to describe the low temperature physics is the boson-boson scattering length, in the present case 
the dimer-dimer scattering length 04. Hence knowing the precise value of this parameter, in terms of the scattering 
length a for fermions making up the dimer, is of utmost importance for the description of this strong coupling limit. 

This problem was first addressed by Haussmann [7] and by Sa de Melo, Randeria and Engelbrecht [8] by different 
methods, which turn out to be equivalent to the Born approximation for this scattering problem. The corresponding 
result is a4 = 2 a. This result was much improved by Pieri and Strinati [5] who considered repeated dimer-dimer 
scattering, instead of a single scattering corresponding to the Born approximation. They proceeded to sum up the 
corresponding series and obtained numerically ~ 0.75 a. Finally an exact treatment was given by Petrov, Salomon 
and Shlyapnikov [101 111] who provided the numerical solution of the corresponding four-body Schrodinger equation. 
This led them 04 ~ 0.60 a. This problem was then taken up by Brodsky, Klaptsov, Kagan, Combescot and Leyronas 
[12] who gave an exact treatment of the same problem, by making use of field theoretical methods. The numerical 
solution of their equations gave naturally the same result 04 ~ 0.60 a as the one obtained by Petrov, Salomon and 
Shlyapnikov. 

In view of the interest in fermionic mixtures made of different elements, such as 6 Li - 40 K mixtures, Petrov, Salomon 
and Shlyapnikov extended their treatment [13) to the case where the fermions making up the dimers have different 
masses. The corresponding extension of the exact field theoretical treatment was provided by Iskin and Sa de Melo 
[14] . who provided numerical results for several mixtures of specific interest. This method has been extended recently 
by Levinsen and Petrov [T5] to the case of narrow Feshbach resonances, aiming specifically at 6 Li - 40 K mixtures. Here, 
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FIG. 1: Diagrammatic representation of the relation Eq.|T| between the full dimer-dimer scattering vertex T4 and the vertex 
$. The black strips are for ladder diagrams corresponding to the dimer propagator T2 given by Eq.|3|. The full line is for the 
t-spin propagator. The dashed line is for the ^-spin propagator. 



in the case of a wide Feshbach resonance, we take over this technique to investigate the large mass ratio analytically. 
Our aim is the same as in our recent work on the fermion-dimer scattering length |16j . namely to gain some insight 
in the dimer-dimer scattering problem which might be used in dealing with more complicated problems where this 
scattering is a building block. 

In the present case we succeed in obtaining such an insight. We indeed find that the Pieri and Strinati [9] approxi- 
mation is asymptotically correct for large mass ratio. Taken with the fact that their result for 04 is also a quite good 
approximation when the two different fermions have equal mass, we come to the conclusion that their approximation 
is quite good for any mass ratio. This is an important simplification in the dimer-dimer scattering problem since this 
means that we do not have to take into account intermediate states where one of the dimer is broken. Except for 
the irreducible process corresponding to the Born contribution, this means that in all the intermediate states in this 
scattering process we have to deal with unbroken dimers. 



II. BASIC EQUATIONS 

As it is quite often done, we will call spin-up and spin-down the two different fermions with respective mass 
and 771^. This is a convenient convention frequently used in cold gases, although there is in principle no link with the 
physical spin of the particles. The dimer is a bound state of one spin-up and one spin-down fermion. Clearly the 
scattering length we are looking for is proportional to the spin-up spin-down scattering length a. For dimensional 
reasons 04 depends only on the mass ratio r = m^/m^. Since exchanging the two fermions does not change the dimer, 
the result is unchanged when r is changed into 1/r. Hence the large mass limit we are mostly interested in can be 
seen equally as the r —> 00 limit or the r — > limit. 

Let us first obtain the generalization for different masses of the equations used by BKKCL [T2]. The two basic 
vertices are again T^{pi,p2] P) and $(qi, q2',P2, P)- They are shown in Fig. [I] The first one describes the scattering 
of two dimers entering with respective four-momenta P+pi and P — p%, and outgoing with four- momenta P+P2 and 
P — P2 ■ In the second one the entering particles are two different fermions with respective four- momenta q\ and (72 
and a dimer with four-momentum P — q\ — q 2 , while the outgoing particles are again two dimers with four- momenta 
P + P2 and P — P2 ■ We will take the convention for $ that the first variable q\ corresponds to the spin- up particle 
while the second one 52 corresponds to the spin-down one, explicitly &(qi t> 92 i-',P2,P)- In contrast with the equal 
mass case $(qi, q^Pi-, P) is no longer symmetric in the exchange of q\ and 172- 

It is clear that the first process in T4 involves the interaction of a fermion of one dimer with a fermion of the other 
dimer. Hence one has first to split open a dimer line into its fermionic components, all the subsequent processes being 
described by This leads to the equality shown in Fig. 1. Algebraically it reads: 

T 4 (p 1 ,p 2 ;P)=J2Gt(k)G i (P+pi-k)$(k,P+p 1 -k;p2,P) (1) 

k 

where Y] = i J dkdu> / (2tt) 4: is for the summation over momentum and energy and the propagators are given by 

k 

Gf^(k) = [lj — fc 2 /2m^ + ie] -1 , with e — > + . In turn we can write an equation for $, in the same spirit as a 
Bethe-Salpeter equation. The first process is a break up of the entering dimer, in order to allow the interaction of 
one of the free fermions with a fermion from the dimer. Actually other interactions between these same fermions 
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may also occur after the first one and the resummation of all these processes leads merely to a dimer propagator T 2 . 
Afterwards any process may follow, which is again depicted by However one must have in mind that by definition, 
in $, the first process can not be an interaction between the two free fermions (otherwise one would merely have an 
entering dimer, instead of two free fermions, which is already accounted for by the entering dimer propagator in T4). 
But such a process is quite allowed after the first interaction, and since they are not accounted for by $ we have to 
add terms to describe them. The repeated interaction of the free fermions leads to another dimer propagator, so at 
this stage we have two dimers. This can be taken as the outgoing state, and the corresponding diagrams correspond 
merely to the Born term for the diagrammatic expansion for $. However it is also possible to have any other process 
after these two dimers have been formed, which is fully described by T4. For each of the process we have described, 
we have actually two possibilities depending on the spin of the involved particles. However this is not true for the last 
T4 term we have just described, since one sees easily that exchanging the spins is equivalent to a change in dummy 
variables. This leads finally to the diagrammatic equation depicted in Fig. [2] which reads algebraically: 

$( qi ,q 2 ;p 2 ,P) = -G t (P - q 2 -p 2 )G i {P - q x + p 2 ) - G X {P - qi - p 2 )G t (P - q 2 +p 2 ) 

-Y J G^k)G i (2P-q 1 -q 2 -k)T 2 {2P-q 2 -k)^(k,q 2 -p 2l P)-Y J G i {k)G^2P- 
k k 

- ]T G t (Q - q 2 )G±(2P -Q- qi )T 2 (2P - Q)T 2 {Q)T 4 {P - Q,p 2 ; P) (2) 
Q 

The dimer propagator is obtained as usual by summing up the ladder diagrams and is given by: 

2"7T 1 

T 2 (P) = = = T 2 (P,E) (3) 

1 ' [i a- 1 - - E) V ' ' W 

for the four-momentum P — {P,E}. Here /i is the reduced mass fj, = ra^m^jM = m^r/(\ + r), while M is the total 
mass M = + = m^(l + r). 

Just as in [12] the scattering length is obtained directly from T4, evaluated at the dimer binding energy Ef, = 
l/(2fia 2 ), the only difference being that we have to use the dimer reduced mass \x and the reduced mass M/2 of two 
dimers: 

2tt \ 2 , , 1% 27r(2a 4 ) 

Since for the scattering length problem we have to take p 2 = and P = {0, — Eb} we do not write anymore explicitly 
these variables. Hence Eq.Q becomes: 

_ _ a 4 _ 7T (1 + r) 5 T 4 (0) 

ai ~ a ~ 2 r-4 (amt) 3 ( j 

In the following we take for convenience a and as unit for length and mass, i.e. we set a = 1 and = 1. This 
last step breaks apparently the invariance of 0,4 under r — > 1/r, but this symmetry is naturally still satisfied in the 
final results. 

A very important simplification is that, just as in |12j . the calculation of the two $ terms in the right-hand side of 
Eq.(2| requires only the knowledge of the "on the shell" value for $ with respect to variable k. This allows to consider 
Eq.(2) only for "on the shell" values for variables q\ and q 2 . We denote <I>(qi,q2) the corresponding "on the shell" 
value of <& (it is not symmetric under the exchange of qi and q 2 ). Similarly when Eq.([2]) is inserted into Eq.([T]), the 
resulting equation requires only the knowledge of $(qi,q 2 ). Finally it is again convenient to continue the equations 
to purely imaginary values for the frequencies. The corresponding value for T^{Q) — T^dQjifl}) = f4(Q,fi) is real 
and the equations contain only real quantities. 

Let us now give the various terms entering the equations for $ and £4. The first two terms in Eq.Q correspond to 
the Born approximation and give to $(qi,q2) a contribution: 

1 + q{ l + qi 

where now q\ = |qi| and = r/(l + r). The two next terms in Eq.Q give two contributions: 

* (a „\- 2 ( rile r2( ~ 1/ ^ - q ^ 2 qi + k ) $ ( qi ' ^ m 

$ $ l(qi,q 2 ) -J^J dk 2/^ + ,f/r + PA + (q 1+q2 +k) 2 (?) 
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FIG. 2: Diagrammatic representation of the equation Eq.|2| for the vertex $. The black strips are for ladder diagrams 
corresponding to the dimer propagator T2 given by Eq.|3|. Full lines are for f-spin propagators. Dashed lines are for ^-spin 
propagators. 



where the integration over the azimuthal angle of k with respect to qi can be explicitly performed, since $(qi,k) as 
well as T 2 depend only on the polar angle between k and qi. This leaves a double integration to be performed. The 
other contribution is 



$*2(qi,q2) 



(2tt) 3 



elk 



T 2 (-l/p - q 2 2 /2r - fc 2 /2, q 2 + k) $(k, q 2 ) 
2/p + qf + ql/r + k 2 + (qi + q 2 + k) 2 /r 



(8) 



where the integration over the azimuthal angle of k with respect to q 2 can again be explicitly performed leaving again 
a double integral. Finally the last term gives, after the change of variable Q — » P — Q: 



$T(qi,q2) 



(2tt)< 



dQ / dQ 



|T 2 (-l/2/i^,Q)| 2 i 4 (Q,r!) 



[l/p - 2tfl + qj + (Q- qi) 2 /r] [1/p + 2^0 + q 2 /r + (Q + q 2 ) 2 ] 



(9) 



Since £4 (Q, VL) as well as T 2 depend only on the modulus of Q, the angular average over Q can be explicitly performed 
by making use of: 



dfli 



1 



1 , ab - fc 2 u.v + A 

= — In 

47T a + k.u b + k.v 2A ab - fc 2 u.v - A 



where 



A = ^{akv) 2 + (bku) 2 - 2abk 2 u.w - fc 4 [(uv) 2 - (u.v) 2 ] 
This leaves again a double integral to be performed. To summarize we have explicitly: 

$(qi,q2) = $s(qi,q 2 ) + $$i(qi,q2) + $*2(qi,q2) + $T(qi,q2) 



(10) 

(ii) 

(12) 
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where the four terms in the right-hand side are defined respectively by Eq.([6|,Eq.([7|,Eq.([8| and Eq.Q. 
Similarly, substituting Eq.Q in Eq.Q we obtain for £4: 

U (q, v) = T B (q, v) + T*(q, v) + T T (q, v) (13) 

After the same change of variable Q — > P — Q the last term reads: 

T T (q, u) = J dClJ dn K(q, v\ Q, Q) \T 2 (-1/2/i - ffl, Q) | 2 U (Q, 0) (14) 

where the symmetric kernel if (q, f; Q, f2) = ^(Q, f2; q, J^) is given by: 

K(q, v- Q, n) = J2 G t(fc)G;(P + 9 - k)G t (k - g - Q)G i {P + Q - k) (15) 

k 

with q = {q, w}. This can be written in a more symmetric way by the change of variable k — > k + (P + q + Q)/2: 

K{q, v; Q, n) = G T (fc + )Gj.( fc)G t (fc + )G i { k) (16) 



After performing in Eq.(15) the frequency integration over ui and making the change k — > k + q, one is left with 



expressions which can be written in terms of products of two rational functions, each one being of the form l/(a+k.u) 



The angular average can be performed by making use of Eq. ( 10 1 . One is left with a simple integration over the modulus 



of k. This means that evaluation of Eq.( 14 ) requires four integrations. However we will not write the lengthy resulting 
expressions, which are only necessary to perform the full numerical calculation for the general case of two different 
masses and mj,. 

Then the Born contribution Tg (q, v) is merely given by: 

T fl (q,i/) = -2Jir(q J i/;0,0) (17) 

Finally the last term T$(q, v) is given by: 

, v 4 f . T 2 (-l//i-fc 2 /2-fc' 2 /2r,k + k') $(k,k') 

T$(q ' V) = ~ J2^ J ^ 2iu + + (k - WM [1/M + 2«/ + k>yr + (k' + q)2] +(^~9) (18) 

Since the result depends clearly only on the modulus of q and the dependence on q is explicit, the integrand can be 
angularly averaged over the direction of q by making use of Eq. ( [To| . This leaves a triple integral over the moduli of 
k and k', and the angle between them. 

In order to solve numerically Eq.(12| and Eq.(13) for <I>(qi,q2) and £4 (q, v) and then obtain 04 from Eq.Q, we 



have discretized these equations. This leads to a set of linear equations, the solution corresponding to a matrix 
inversion. This has been performed making use of the standard LAPACK routines. This works quite well for not too 
high mass ratios, and we estimate the precision of our results to be typically 2%. However for mass ratio typically 
above 500, the numerics becomes unreliable. This is easily understood from the somewhat singular features which 
emerge from our anaytical solution, presented below, for very large mass ratio. We note finally that we never find 
any zero eigenvalue for the matrix to be inverted. This means that the bound states discussed in |13j do not play any 
role in the calculation of 04. 

III. SUM RULE 

In a way completely analogous to what we have found in our study of the atom-dimer scattering length [TB] , a quite 
useful sum rule can be obtained by analyzing the convergence of the various integrals when the variables go to infinity. 
It is natural to assume that the solutions $(qi, 02) and t± (q, v) we are looking for have physical ranges corresponding 
to their variables, and that beyond these ranges these functions go rapidly enough to zero for the various integrals 
to be convergent. This assumption is confirmed by our numerical calculations. This property allows to study the 
behaviour of $(qi, q2) and t± (q, v) for large values of the variables. 



Let us first consider the <4 equation Eq.(13) and begin by the Born contribution Ts(q,^) given by Eq.(17) and 
Eq.(15|. Since we are interested in the case where q and v are large, P is negligible in this regime and we are left 
with: 

T B (q,u) « -2j2Gt(k)G l (q-k)G t (k-q)G^(-k) (19) 



() 



It is then convenient to consider more appropriate energy variables. Since energy is homogeneous to momentum 
squared, we set v = q 2 and similarly for the integration variable uj = For large values of q and v, the natural 
variable to consider is p = \/q 2 + q 2 ., corresponding to introduce a radial coordinate for |q| and q v . Similarly we can 
introduce r — -\/k 2 + /c 2 (not to be confused with the mass ratio). Now for example G^(k) is homogeneous to *■ J 



and similarly for the other Green's functions entering Eq. ( 19 1 . Hence for homogeneity reasons we have Tg (q, v) ~ p 3 
since the summation J"\, introduces a factor homogeneous to r 5 . More precisely, making the change of variable r = fp 



allows to obtain in Eq.(19) a prefactor p~ 3 , the remaining factor being independent of p. One can naturally write 
more explicit expressions for Tg(q, v), although it is complicated to go to a full analytical result. However we will not 
need these expressions and they merely confirm our result that Te(q, v) ~ (q 2 + q 2 )~ 3 ^ 2 which comes simply from 
our homogeneity analysis. 



We can now make a similar analysis for the term T^(q, v) in the £4 equation. Indeed since in Eq.(14) we are 
interested in large values of q and v, while the integration variables Q and Q are effectively bounded by the factor 
ti(Q, £1), the kernel reduces in this limit to: 

K(q,u;Q,n)^K{q,p;O,0) (20) 



as it is quite clear from the symmetric form Eq.(16l for K. Hence this kernel factorizes out and we are left with: 

r r (q, v) = vAj^Cq, f ; o, o) J dQ J dn \t 2 (-1/2/, - in, q) | 2 1 4 (q, n) (21) 

which has naturally the same behaviour as the Born term, namely Tr(q, v) ~ (q 2 + q 2 )^ 3 ^ 2 = p~ 3 - 

Finally we analyze in the same way the last term T$(q, v) in this t± equation from its explicit expression Eq.(18|. 
With the k and k' variables bounded by the $(k, k') factor, the q and v dependence comes explicitly from the 
denominators and is given by I#(q, v) ~ (q 2 + q 2 )~ 2 = p~~ A . Hence it converges toward zero faster than the two other 
terms and the overall behaviour of t<i(q, v) is apparently ^(q, v) ~ p -3 . 



However when we insert this behaviour in the integral factor found in Eq.(21 1 for the large (q + q 2 ) behaviour we 
find a discrepancy. Indeed from Eq.|3|, and for large Q and n, we have \T 2 (—1/2// — in, Q) | 2 ~ R~ 2 , where we have 
again made the change = for the energy variable and introduced the radial variable R = \/ Q 2 + Q 2 ? . This leads 
to |T 2 (— 1/2/x — ifl, Q) | 2 t 4 (Q, f2) ~ R~ 5 . However we have dQ dfl ~ R 4 dR so that the integral in Eq.(21 ) diverges 
as dR/ R. This is in contradiction with our initial assumption that £4 (Q, il) insures the convergence of the integral. 

The only escape is that the contributions from the Born term and from the Tt, which have exactly the same power 
law dependence, cancel out so that the actual decrease of (Q, il) is faster than R~ 3 . This occurs if the corresponding 
coefficients cancel exactly, which leads to the sum rule: 

j^- J dctj dn \t 2 (-i/2 M - in, q) 1 2 u (Q, n) = 2 (22) 

Upon checking this result numerically, we have found that it is in very good agreement with our calculations. 

Let us now analyze in the same way the $ equation Eq.(12|. First the Born term is explicit and, for large qi and 
q 2 , it behaves as: 

*s(qi,q 2 )«-8// 2 ^ (23) 

9192 

Then we consider the 3>T(qi, q 2 ) term given by Eq.([9|. Again the factor t± (Q, Q), assumed to go rapidly enough to 
zero for large Q and f2, makes these variables Q and f2 effectively bounded. Hence when we consider very large values 
of qi and q 2 , we can forget about Q and Q in the explicit denominators, and the product of these denominators is 
given in this case by q 2 (l + l/r)q 2 (l/r + 1) = ((/i^/m) 2 - This leads for large q\ and q 2 to: 

4 u 2 f f°° 

Mqi,«fc)« 7^4i / rf Q / dn\T 2 (-i/2p-in,q)\ 2 u(Q,n) (24) 

(2tt) 4 qfq^ J 

We consider finally the $$i and $$ 2 terms. Our assumption is that $(qi,q 2 ) goes rapidly enough to zero at large 
qi and q 2 for having the integrals in $$i and <I>$ 2 converge for large k due to the $ factor. This means explicitly 
that J rfk$(qi,k) and Jdk$(k, q 2 ) are convergent. However this is not true if we take the behaviour ~ 1 / (gf Q§) 
produced by the Born term and the $t term, since this gives an integral J dk/k 2 ~ J dk, which is divergent. If the 
dominant behaviour was coming from the terms themselves, then $(qi,q 2 ) would go even more slowly to zero at 
infinity and the situation would be even worst (one can actually check that this case does not arise). Hence we have 
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again a contradiction with our starting hypothesis, but there is also the same way out, namely that the contribution 
of the Born term and the $t cancel precisely, which leads to a decrease of $(qi, q2) faster than ~ l/C^iQ'I ) f° r large 



qi and c\ 2 - We see from Eq.(23| and Eq.(24) that the condition for this to happen is exactly the sum rule already 
found Eq.(22|. This result is not so surprising when we remember that the £4 equation is obtained by carrying the $ 
equation into Eq.Q. 

We note finally that, in our work on the atom-dimer scattering length |16j . we have shown that the sum rule we have 
found in this case is merely a direct consequence of the Pauli principle, namely the fact that two identical fermions can 
not be found at the same place. In the present case such a physical interpretation is less obvious since the frequency 
dependence of £4 (Q, f2) enters. Anyway we have not tried to find a specific physical interpretation. 

IV. VERY HEAVY MASS EQUATIONS 

Let us now consider how the preceding equations simplify when we consider the limiting case where the two masses 
and are very different. We can equivalcntly consider that the mass ratio r = m^/m^ goes to zero or to infinity. 
We will take this last option in the following since it is somewhat more convenient with the asymmetric notations we 
have chosen. This implies fi — >■ 1 in this limit and the Born term Eq. ([6]) in the $ equation becomes: 



*s(qi,q 2 ) = -8— j — 5 (25) 



1 1 

1 + qf YTq'i 

In the $$1 and <E>$2 terms we may replace in Eq.([3]) T 2 (E,P) by 2w/(l — \J — 2E) since M — > 00 in this limit. This 
leads to: 

11/" cf>°° (rt ]r\ 

- -^tttw^ + (* + 1 + (26) 



and 



If 1 <E>°°(k a 9) 

•»(*■*)— 5?/* vm^i 5+3^ 



Let us finally consider the $t term which gives rise to more problems. If we proceed in the same way as above, we 
find: 

4 f f°° 1 t 00 (O ft) 

m«.«)- V g ji f* i L da FmiFwra wrmdip (28) 

However we see that the integration diverges as J dtl/fl 2 for O — >• 0, since + 2ifl — 1| 2 ~ 4J1 2 in this limit. 
Indeed the explicit denominators go to finite values and there is no reason to have ti(Q, f2 = 0) = 0. This can anyway 
be checked numerically. In particular this does not happen for Q = since, in the large r limit we have from Eq.([5]) 
(with our reduced variables): 

a 4 = |rt 4 (0,0) (29) 

The existence of this divergence shows that we have handled the \T 2 \ 2 term too rapidly, since it is at the origin of 
the divergence which does not exist naturally in the general equation Eq.([9]). We must keep M w r without setting 
immediately l/M — 0. In this case we have to write: 

9tt 

T 2 (E, P) = ; (30) 

V ' 1 - y/-p 2 /r-2E 



which leads, instead of Eq.(28), to 



4 f f°° tf (Q,n) |VmgT+QVr + lj 

(2tt) 2 J [l-2^^2 + g 2 ][l + 2^^! + (Q + q 2 ) 2 ] 4fi 2 + Q 4 /r 2 



We see that the divergence has disappeared. Instead we find a factor 1/(4S! 2 + Q 4 /r 2 ) which, in the limit r —> 00, is 
strongly peaked around 57 = and acts in practice as a 5 function: 

I 7T7* 

' -6(Sl) (32) 



4fi 2 + Q 4 /r 2 2Q 2 
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As a result only ti(Q, Q = 0) appears in the $ equation so that we may write the equation for £4 only for fl = 0, 
which is a quite convenient simplification. Taking into account that i4(Q,f2) depends only on the modulus Q = |Q|, 
we set: 



U{@) =rt 4 (Q,O = 0) 
which is just the quantity coming in the scattering length we are looking for: 



0,4 



h{0) 



(33) 



(34) 



Taking the limit r — > 00 in the numerator of Eq.(31) and performing the angular average over the direction of Q, wc 
finally end up with: 



$t (qi.qa) = 



(l + «l)<72 Jo 



dg ^M ln i + (o + ?2 ) 2 



Q 



1 + (Q - g 2 ) 2 



(35) 



We can now write the £4 equation Eq.(13l with the variable v set to zero, which is the only thing we need as we 



have just seen. The Born contribution is given by: 

Ts(q,0) 



-2if(q,0;0,Q) 



(36) 



where K (q, 0; 0, 0) is obtained from Eq.( 15 ). Since in the limit — > 00, reduces to [w + ie] 1 the calculation 

is fairly simple and leads to: 



1 



7T q 2 + 4 



(37) 



The 7$ term is also easily obtained in the limit r — > 00 from Eq.(18), making use of the limiting expression for T 2 

l + {k' + q) 2 



and performing the angular average over the direction of q. This gives: 

k 2 1 



T|°(q,0)^^( g ) = ^ 1 

^ q Jo 



dk 



1 + k 2 y/2T¥ - 1 



dk k In 



i + (fc'-<z) 2 



dn^ 

47T 



$°°(k, k') 



(38) 



Finally in the last term of the <4 equation, coming from Eq.(14|, we meet the same troubles as in the last term of 
the $ equation if we use the simple limiting expression for T2. In the same way as what we have done to derive 
(Qi> Q2), we have to use Eq.(30) for T2. This leads in the same way to a factor proportional to S(fl) justifying 



the fact that we write the t 4 equation only for zero frequency. Actually the last two factors in Eq.pT| are the same 
as those appearing in Eq.([9|. Hence the only difference is that we have now to calculate the kernel K(q : 0; Q, 0) for 
r — > 00. This proceeds just as for the Born contribution Eq.(37) and leads basically to the same result, provided q is 
replaced by q + Q. This gives: 



tf(q,0;Q,0) = 



1 



7T (q + Q) 2 + 4 

When the angular average over the direction of q is performed, just as in the preceding term, one finds: 



23?(q,0)=t§?(g) 



(Q + q) 2 



«~qJo ^ V ln A + (Q- q y 
To summarize the $ equation becomes: 

<f°° ( Ql , q 2 ) = $~ ( qi , q 2 ) + $~ ( qi , q 2 ) + ( qi , q 2 ) + $°? ( qi , q 2 ) 



(39) 



(40) 



(41) 



where the four terms in the right-hand side are defined respectively by Eq.|25j), Eq.(26), Eq.(27) and Eq.(35), while 
the t\ equation reads in this limit, with the definition Eq.(33): 



-3° («)=*£ («)+«§? (?) 



(42) 



where the three terms in the right-hand side are defined respectively by Eq. ( 37 1 , Eq. ( 38 ) and Eq. ( 40 1 . 

Although these equations correspond to a very important simplification with respect to the original ones, they are 
still too complicated to be solved analytically as such. In the following section we will show that they can be further 




FIG. 3: (Color online) Dimer-dimer scattering length 04 as a function of the mass ratio r (logarithmic scale). The red full line 
is the exact numerical result obtained from Eq.(12l and Eq.(13l. The green dotted line is the numerical result obtained from 
Eq.pUl and Eq.p2j). The blue dashed line is the result obtained from Eq.(|42| by setting $°°(qi,q 2 ) = 0. 



simplified, leading to an analytical answer. However it is of interest to solve them numerically to obtain the scattering 
length a 4 and to compare the result to the exact numerical solution of the original equations Eq.(12) and Eq.(13l 



This is done in Fig. [3j It is quite surprising to see that, already for a mass ratio slightly above 10, the result from 
these asymptotic equations coincide with the exact one within numerical precision. This provides naturally a further 
validation of these asymptotic equations. 

Finally let us consider the sum rule Eq.(22) in this limit. We have to handle the 
done just above. This implies in the same way that only i 4 (Q, £1 = 0) 
with the very simple relation: 



T 2 | 2 factor carefully, as we have 
appears in the sum rule and we end up 



dQ*4°(Q) = 1 



(43) 



V. DISCUSSION OF THE VERY HEAVY MASS LIMIT 



Although the equations have been much simplified in this heavy mass limit compared to the general ones, they 
are still fairly complicated. We will nevertheless be able to come to a very simple conclusion by showing that a 
quite natural hypothesis on the behaviour of if(Q) is fully consistent with the equations, and is in agreement with 
results found numerically. However being able to prove that this is the only possible solution looks a very difficult 
mathematical problem. 

The natural hypothesis stems from the sum rule Eq. ( 43 ) and from the fact that the scattering length a 4 grows when 
the mass ratio increases, as it is known from preceding work [111 I14| and from our own numerical calculations. It is 
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natural to assume that 04 grows indefinitely. From Eq.(34| this means that ^(O) grows indefinitely for large mass 
ratio. However the sum rule Eq.(43) puts a constraint. If we assume, as we have already done, that if(Q) decreases 



rapidly when Q is large and has a fairly regular behaviour, the increase of ^(0), with fixed surface under the curve 
if(Q) forced by the sum rule, implies that becomes very peaked around the origin for very large mass ratio. 



In such a case we can further simplify the equations. Let us take first take <&5?(qi,q2) given by Eq.(35). Since 



tf(Q) is peaked around the origin, Q is actually forced to be small. We can then expand the logarithm into 
ln[l + (Q + q 2 ) 2 ]/[l + (Q- q 2 ) 2 } =s 4Qq 2 /(1 + q\). This leads to: 



+ + Jo 



dQtf{Q) = 



(l + q 2 )(l + q 2 ) 



from the sum rule Eq.(43). Hence 3>2?(qi,q 2 



cancels exactly the Born contribution Eq.([25|) 

>Q2 



(44) 



Therefore the only 



contributions left in the right-hand side of the $ equation are $^i(qi, q 2 ) and ^^(qi, 12)- However this means that 
the $ equation is now an homogeneous linear equation in <5>°°(qi, q 2 ), without source term. Barring a singular kernel 
(which would make impossible in general to solve this $ equation), the only solution is merely $°°(qi, q 2 ) =0. 
Hence we are only left with the £4 equation Eq.(42) from which the t^(q) term has disappeared. In the equation 



Eq.(40l we can expand, for the same reason as above, the logarithm into ln[4+ (Q + q) ]/[4+ {Q — q) ] — AQq/(A + q 2 ) 



which gives: 



7T 4 - 



dQt?(Q) 



7T 4 - 



(45) 



again from the sum rule. And we see again that this term t^(q) cancels exactly the Born term Eq.(37). Hence the 
£4 equation Eq.(42) becomes merely tf(q)/r = 0. But this is perfectly consistent with the fact that we deal with 
r = 00 limit. In conclusion we find that the equations are perfectly satisfied in this limit by a function t£°(Q) very 



strongly peaked around the origin and satisfying the sum rule Eq.(43). The only trouble is naturally that we have 



been unable to extract any information. In order to obtain results we have to take more carefully advantage of the 
fact that if{Q) is strongly peaked. 

For this purpose let us go back to the $ equation Eq.(41). The source term (qi,q 2 ) + $jf(qi,q 2 ) can be 
factorized into: 



*B(qi,q3) + $T(qi>qa) 



(1 + 9?) 



S(q2 



where: 



S{q2) = — dQ -1 , 7s~i y) 2 



(46) 



(47) 



In the above discussion, valid for an extremely peaked function if(Q) ~ S(Q), we had S , (g 2 ) = 0. When r is large, 
but finite, if{Q) is peaked around the origin which makes Q effectively bounded. When g 2 is large, we can again 
expand the logarithm in Eq.(47) as above and reach again the conclusion that 5(g 2 ) = 0. Hence S , (g 2 ) is also peaked 
around the origin. 



Moreover, even if S(q2) is not zero, it retains the following exact property: 

dkk 2 S(k) = 



Indeed making use of the sum rule Eq.(43) we can write from Eq.(47): 

k , l + (Q + k) 



dkk z S{k) 



dQtf{Q) / dk 
Jo 



4k 2 



Q h \ + (Q- fc) 2 1 + fc 2 



(48) 



(49) 



The integral over k can be calculated analytically, and it is found to be zero whatever the value of Q. 

These two properties imply that the equation has, to a very good approximation, a factorized solution of the 
form: 



$°°(qi,q 2 ) = F( qi )S(q 2 ) 



(50) 



Indeed when this factorized expression is substituted in the $ equation Eq.(41 ) we see that we can factorize S(q 2 ) not 



only in the source term Eq.(46), but also in 3>|? 2 (qi,q 2 ) (see Eq.(27)). In the remaining term $$i(qi,q 2 ) given by 
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Eq. p6| , when we substitute Eq. ( 50 1 , the factor S(k) prohibits large values for |k| , as we have just seen. But when |k| is 
small, or at most of order unity, it is a very good approximation to neglect it in the denominator 2 + q 2 + (qi +q2 +k) 



Hence we are left in this term with the integral J dk.S(k) — An dkk 2 S(k) — as we have just seen in Eq.(48|. 
Therefore the contribution of the $^ term is zero to a very good approximation, which shows that the solution of the 
$ equation is indeed of the form Eq.(50). We have checked that this decoupling is properly satisfied asymptotically 



by the results of our numerical solution. 
The resulting equation for F(q) is: 

F(q) 



1 



dk 



F(k) 



+ W - 1 2 + k 2 + q 2 



(51) 



This equation is easily solved and gives a well behaved solution. This is most easily understood if one notes that 
Eq.(51| has an interesting physical interpretation. Indeed, as we have seen, $(qi,q2) describes the scattering of two 
free fermions on a dimer. In the large mass ratio limit, it is reasonable to assume that the free heavy fermion does 
not play any role and that only the scattering of the light fermion on the dimer is relevant. In this case we are back 
to a problem first investigated by Skorniakov and Ter-Martirosian |17j for the case of equal masses, and for which we 
have recently found an analytical solution in the case of very different masses [TB] . However there is a slight difference 
between our case and the fermion-dimer scattering length problem. In this last one, since the kinetic energy is zero, 
the total energy is just the dimer binding energy — E^. In our case, since by definition of $ the final state is made of 
two dimers with zero kinetic energy, the total energy is —2Ejj, that is twice the dimer binding energy. The two terms 
2 present in the integral in the right-hand side of Eq.(51 ) can be tracked back to this total energy, taking into account 



that we have used reduced units. If we really had a fermion-dimer problem, these two 2 should be replaced by two 1. 
In this case, making this substitution, we would rather have the equation: 



F'(q) 



2 f°° 

n Jo VI 



F'(k) 



k 2 -l 1 + k 2 



Making the change F'(q) — 2a^(q) / 1 + q 2 + 1) we end up with: 



a3{q) 



dk 



a 3 (k) 



1 + k 2 



(52) 



(53) 



which is exactly the equation we had obtained |16j when one fermion in the dimer is very heavy while the two other 
ones are very light. We have found the analytical solution a 3 (q) = 1/(1 + q 2 ). This makes it easy to understand that 
the solution of Eq.(51) is very similar and indeed we have found the analytical solution: 



F(q) = 



(i + « 2 )(v^ 



1) 



(54) 



This analysis makes it also possible to understand physically the factorization Eq.(|50[). Indeed since the evolution of 
the heavy and of the light fermions decouples, it is quite natural that in the vertex $, which describes their evolution 
in the presence of a dimer, their contributions factorize as it is the case for the wavefunction of two independent 
systems. 

Having seen that the solution Eq.([50| for $ is perfectly acceptable, we can substitute it in the equation for 
Eq.(42|. However in the term t^(q), given by Eq.(38), we can make use of the fact that we are only interested 



in small values of the variable q since we know that t^°(q) is peake d a round the origin. In this case writing again 
ln[l + (k' + q) 2 ]/[l + {k 1 - q) 2 } ~ 4k'q/(l + q 2 ) and substituting Eq.([50} leads to: 



**(?) 



1 



1 



dk 



F(k) 



i + fc 2 VzTk 2 



poo 

- / dk'k' 2 S{k') = 
1 Jo 



(55) 



where the last equality comes again from the property Eq. ( 48 ) . Hence <!> disappears entirely from the equation for t± 
and we are left explicitly with: 



1 



if(q) 



4 
7r q 



1 

2 j_ 4 



1 1 

7r q 



dQ^ 



'(Q) , 4 + (Q + g) 2 
Q n A + (Q-q) 2 



(56) 



We have checked numerically that taking $°°(qi,q2) = in Eq.(41| and Eq.(42) leads to the correct solution. The 
exact numerical solution of Eq.(42) with <&°°(qi,q2) = leads for to the results displayed in Fig. [3j We see that. 
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for a mass ratio above 10, it gives a fairly good result for 04, and that when increasing further the mass ratio the 
result converges toward the exact one. We note that the disappearance of $ from the equation for t± means that 
the approximation made by Pieri and Strinati [S] in their work on the dimer-dimer scattering length is fully valid in 
the limit of large mass ratio. Indeed they made a ladder approximation, neglecting the possible breaking of a dimer 
described by $ and retaining only repeated scattering between the two dimers. We see that, since for the equal mass 
case, this approximation gives 04 ~ 0.78 instead of the exact [THIUS] 04 ~ 0.60 and since it becomes exact in the limit 
of very different masses, it turns out to be quite a good approximation for any value of the mass ratio. 
If we perform as a final step the change of function tf(q) = irif (q), Eq.(56) becomes: 



r 



(?) 



1 



1 



1 1 

2tt q 



dQ 



tf{Q) , l + (Q + qY 



Q 



In 



1 + (Q - qf 



(57) 



This equation is identical to the one we had in the fcrmion-dimer problem [TB] when we considered the large mass 
ratio domain. The only apparent difference is that, in the left-hand side of the equation, the coefficient in |16| is 
m^/m^, while in the present case it is 1/r = m^/m^. However since in |16) . the considered limit is m^/m^ — > while 
we consider here m^/m^ — > 00, and since the scattering length 04 is unchanged under the exchange m-j- f-> toj., the 
equations are indeed identical. 

Hence we have tf (0) = 0,3 — /a in terms of the result obtained in [TB] . Using Eq. ( 34 ) and making the exchange 
rn^ O mi, we obtain our final result for the value of the dimer-dimer scattering length in this large mass ratio limit 



03 

(24 = — - = 



ln(m t /m;) - ln(ln(m t /m;)) + 2C 



(58) 



where C — 0.577215... is the Euler constant. On the other hand for m-j- = 771^ we have 04 ~ 0.60 while 03 ~ 1.18. 
Hence in this case we also have 04 ~ 03/2. Hence it can be guessed that the equality 04 = 03/2 is approximatively 
valid whatever the mass ratio. Indeed the two quantities 04 and a$/2 are displayed in Fig. [4] and it can be seen 
that their difference is at most ~ 0.1, obtained for ln(ra^/m^) ~ 2. This difference is naturally expected since a 4 is 
invariant under O while this is not the case for 03. 

The last expression in our result Eq.(58), which is displayed in Fig. |4j is in agreement with the work of Marcelis, 
Kokkelmans, Shlyapnikov and Petrov (MKSP) [18] for 04 in this large m^/m^ regime. They addressed this problem 
with the 4-body Schrodinger equation, which they solved in this regime by a Born-Oppenheimer approximation. They 
found the approximate relation 2S4 — 2C + ln(m^/2m^a4). To dominant order (i.e. omitting the 2<i4 in the right-hand 
side) it gives 2S4 = 2C + ln(TO|/77ij.) in agreement with two terms of our Eq.( |58[ ). Our third term is recovered by 
inserting this expression in the right-hand side of their relation, corresponding to the next step in a recursive solution 
of this equation, and keeping only the dominant contribution. 



In this spirit one can find an analytical formula, slightly different from Eq.(58l but equivalent for large m^/rn^, 
which is much closer to the exact numerical result than our Eq. ( 58 1 or also than the exact numerical solution (also 
displayed in Fig.|3| of MKSP approximate relation. In Eq.(58) the trouble comes at low m^/mj. from the divergence 
of ln(ln(m-t-/m^)) when m^/m^ — > 1. This is cured by introducing a constant in the logarithm. This constant could 
be adjusted for fine-tuning. But taking it equal to 1 is both simple and in very good agreement with numerics. This 
leads to: 



04 



ln(m t /m;) + 2C - ln[ln(m t /m^) + 1] 



(59) 



As it can seen from Fig. [4j it is fairly close to the exact numerical result for 04 (and even extremely close to 03/2). 
In general the agreement is expected to be limited by the fact that the exact result has the property to be invariant 
under m-f O TO4. (which leads to a zero derivative with respect to ln(m^/m^) for m^/m^ = 1), while an approximate 
analytical expression will not have this property. However we see that Eq.(59) has precisely a zero derivative with 
respect to ln(m-j-/m^)) for m^jm^ = 1. 



VI. CONCLUSION 



In this paper we have studied the dimer-dimer scattering length 04 for a two-component Fermi mixture with different 
fermion masses and respectively. For this purpose we have made use of the exact field theoretical method 
already present in the literature |12l 1141 115] . The corresponding equations have been solved numerically for any value 
of the mass ratio m^/m^. However our main aim has been to study the large mass ratio domain. In this range we have 
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FIG. 4: (Color online) Dimer-dimer scattering length 04 as a function of the mass ratio r (logarithmic scale). The red full line 
is the exact numerical result obtained from Eq.(12l and Eq.(13|. The purple dotted-dashed line is the exact numerical result 
for 03/2. The black dashed line is Eq. 
numerical solution of the equation 2<i4 



58|| The blue dotted line is Eq.([59j). The light blue double-dotted-dashed line is the 
2C + ln(mt 727714.04) found by MKSP [18]. 



been able to simplify the equations enough to obtain an analytical solution. More specifically we have shown that 
our final equation is essentially the same as the one obtained in the fermion-dimer scattering problem with scattering 
length 03. In this way we have shown that, for large mass ratio, we have the very simple result 04 — Since this 

relation is also correct with a very good precision for = m^, ~ 03/2 happens to be valid for any mass ratio with 
a quite reasonable precision. We have found for the scattering length an analytical formula which implements this 
finding and gives the result with a very good precision for any mass ratio. For the large mass ratio domain, our result 
agrees for the dominant terms with the earlier one [18] obtained from a study of the 4-body Schrodinger equation. 

A very important finding in our investigation is that, in the large mass ratio domain, the dominant process in the 
dimer-dimer scattering is the Born contribution, with resummation of any number of such processes. Other more 
complex processes, where for example one fermion propagates freely while the other one scatters on the other dimcr, 
become negligible. Retaining only these repeated dimer-dimer Born scatterings is precisely the approximation made 
by Pieri and Strinati [5] in their study of the dimer-dimer scattering. More precisely they made this approximate 
treatment in the case where the fermion masses are equal m-f = m\.. In this case there is no general justification for 
this approximation. Nevertheless it gives for the scattering length a result 04 ~ 0.78 a which is not so far from the 
exact one 04 ~ 0.60 a (mostly if one keeps in mind the simple Born result — 2 a). As a consequence we come to the 
important conclusion that the Pieri and Strinati approximation is a fairly good one whatever the mass ratio. This is 
quite interesting since the processes they retain are much simpler than the ones which have to be considered in full 
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generality. 
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